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Abstract 

The inherent structure approach, wherein thermodynamic and structural 
changes in glass forming liquids are analyzed in terms of local potential en- 
ergy minima that the liquid samples, has recently been applied extensively 
to the study of thermodynamic aspects of glass forming liquids. The evala- 
tion of the configurational entropy, which arises from the multiplicity of local 
energy minima, plays a central role in such analyses. Results are presented 
here concerning the calculation of configurational entropy based on computer 
simulations of a model liquid, and the application of the inherent structure 
formalism to the study of the glass transition locus, and the fragility of glass 
forming liquids. 

I. INTRODUCTION 

Most liquids, when cooled to low temperatures in such a fashion that the equilibrium 
crystal phase does not form, undergo a lower temperature transformation to a solid phase 
which is amorphous, i.e., lacking the long range periodic structure that characterizes crys- 
talline solids. The amorphous phase is a glass, and the transformation process is referred to 
as the glass transition. In the process of cooling towards the glassy state, liquids display a 
rapid increase in viscosity, and display other unusual dynamical features such as stretched 
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exponential relaxation and heterogeneous dynamics |T]-§[. The glass transition observed in 
experiments (the "laboratory glass transition") is understood to be not a transition in the 
thermodynamic sense but a falling out of equilibrium of the liquid. Whether a thermo- 
dynamic transition underlies the laboratory transition, and how an explanation of such a 
thermodynamic transition is linked (or not) to a detailed understanding of the unusual dy- 
namics displayed by liquids at low temperatures, is a topic is considerable current research 
||. A particular approach that has been studied, based on the analysis of local minima of 
the liquid's potential energy (termed "inherent structures") || is described here. Section 
H outlines the approach based the study of inherent structures. Section |TTT] describes the 
details of computer simulations of a model liquid employed in implementing the analysis 
of some properties of glass forming liquids, along with the results. Section [TV] contains a 
summary of the results and conclusions. 



II. THE INHERENT STRUCTURE APPROACH 

The disordered structure of a liquid has the implication that the energies of interaction 
between particles will generally be very complicated, and that the part of configuration 
space explored by the material in the liquid state is characterized by the presence of many 
local minima of the potential energy. Such is the case also, e. g., in a crystalline solid if one 
allows for the presence of defects. The use of the phrase "energy landscape" to describe the 
complicated interactions in a glass forming liquid (and other disordered systems) therefore 
contains in addition the expectation that the complicated potential energy topogrphy plays 
an essential role in determining the properties of the system. If such is the case, it is desirable 
to attempt a description of glass forming liquids in terms of quantities that define the nature 
of the potential energy landscape. In the inherent structure approach, one considers the 
decomposition of the 3N dimensional (for an atomic liquid) configuration space of the liquid 
into basins of individual local potential energy minima. A basin of a given minimum is 
defined as the set of points in the configuration space (or configurations) which map to that 
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minimum under a local energy minimization. The canonical partition function of the liquid 
can then be expressed as a sum over inherent structure basins, the summand being partial 
partition functions defined for individual basins. In the following, the equations are written 
for a two-component atomic liquid, since the model liquid described in the next section is 
such a liquid. In turn, the sum over basins is written in terms of (a) a distribution of minima 
in energy, and (b) the free energies of basins, as follows: 

Q(N, p, T) = A~™ ^-L^ J dr N exp [-/*&] (1) 

= ]T exp [-/3$ a ] A- 3N [ dr N exp - 

= j d<$> exp[-P(<Z> + F vib (<f>,T))} 

= Jd<$> exp [-/?($ + F vib (<S>, T) - TS C ($))} 

where <£> is the total potential energy of the system, a indexes individual inherent structures, 
$ a is the potential energy at the minimum, f2($) is the number density of inherent structures 
with energy $, and the configurational entropy density S c = ks InQ. The basin free energy 
Fvib($a, T) is obtained by a restricted partition function sum over a given inherent structure 
basin, V a . A is the de Broglie wavelength, Na and N B are the number of A and B type 
atoms in the two component liquid, T is the temperature, and p the density of the liquid. 
In the following, the dependence on p is not explicitly stated always since the interest is in 
T dependent behavior at constant density. The configurational entropy of the liquid arises 
from the multiplicity of local potential energy minima sampled by the liquid at temperature 
T, and is related to the configurational entropy density above by 

S C (T) = |rf$5 c ($)P($,T), (2) 

where 

P($, T) = exp + F vib ($, T))] /Q(N, p, T), (3) 

= exp + F vib ($, T) - TS C ($))} /Q(N, p, T), 

is the probability density that inherent structures of energy $ are sampled at temperature 
T. In the above expression for the partition function, an assumption has been made that the 
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basin free energy does not differ for different basins of the same inherent structure energy. 
Without reference to the distribution of minima, the configurational entropy can be defined 
as the difference of the total entropy of the liquid and the vibrational entropy of typical 
minima sampled at a given temperature: 

S c (p, T) = S tota i(p, T) - S vib (p, T). (4) 

The "entropy theory" of Gibbs, Di Marzio and Adam |7]|| define the ideal glass transi- 
tion, underlying the laboratory transition, as an "entropy vanishing" transition where the 
configurational entropy vanishes (the configurational entropy is not, however, defined in pre- 
cisely the same way in |7]||). A similar picture also emerges from the study of mean field 
spin glass models and calculations motivated by them pHTTf. Whether such a transition 



exists for real materials is still a matter of debate |12| , [13| . The calculations below produce 
such an entropy vanishing transition but it must be kept in mind that they result from 
extrapolations which may not be valid. 

Further, Adam and Gibbs || theory relates the configurational entropy to relaxation 
times in the liquid: 

A 



t = r exp 



(5) 



where A is a material specific constant. The validity of this relation has been verified by 
numerous experimental studies (which typically use the excess entropy of the liquid over 
the crystal in place of S c ) and computer simulation studies fi~4]-|i~7| (where configurational 
entropy is evaluated). Further, if S c has the form TS C = Kag(T/Tk — 1), the Adam-Gibbs 
relation results in the Vogel-Fulcher-Tammann-Hesse (VFT) form, observed to describe the 
T dependence of viscosity, as well as diffusivity and relaxation times, in many glass formers. 
The VFT relation may be written as 

1 



t = To exp 



(6) 



K VFT (T/T - 1 

where T is the temperature of apparent divergence of viscosity. Kvft is a material specific 



parameter quantifying the kinetic fragility. Fragility is a measure of how rapidly the viscosity, 
relaxation times, etc of a liquid changes as the glassy state is approached ||19|| . 

That the basin free energy F v n, arises from "vibrational" motion within individual basins 
is emphasized by the suffix vib. If this motion is sufficiently localized around the minima, a 
suitable procedure would be to approximate the basins as harmonic wells, and to evaluate 
the basin free energy within this approximation. The validity of such a procedure has been 
tested recently in various studies |^J^JlO|JT6|JT7|j2^] . It is found that below the temperature 



where the liquid begins to exhibit aspects of slow dynamics, (non-Arrhenius behaviour of 
relaxation times, and stretched exponential relaxation) p3| , p0| , p2| , a harmonic approximation 
of the basins is reasonable. Additional evidence is provided in this paper concerning the 
validity of this approximation, as well as a classical (as opposed to quantum mechanical) 
treatment of the basin partition function. A classical calculation of the basin free energy 
yields 

3N hv- 

F mb = k B Tj2^-^-, (7) 

i=l k B± 

or equivalently, the basin entropy, 

(8) 

where are the vibrational frequencies of the given basin, and h is Plank's constant. From 
the form of S v n, it is apparent that the entropy difference between two basins arises solely due 
to the difference in their frequencies. Thus, such entropy differences remain finite as T — > 
which is unphysical as the basin entropy of each basin and therefore their difference must 
go to zero for T = 0. The extent to which this artifact affects the applicability of classical 
calculations employed in many recent studies depends on the deviation of the quantum 
mechanical result at the glass transition temperature (since below the glass transition, the 
system occupies a single basin). The quantum mechanical basin free energy is given, in 
terms of the same vibrational frequencies z/j, as p4| 



3N 

F vib = k B T hi 
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exp — — — — exp 
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(9) 



and the basin entropy as 
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A comparison will be made in the next section of these two (Eq.(|8|) and Eq.([TT])) expres- 
sions for S V ib- Calculations based on Eq.@, where the vibrational frequencies are obtained 
numerically for energy minima generated in simulations, indicate |l7j (see also [^,^|) that 
the difference in S V i b , between basins is roughly linear in the basin energy Thus one can 
write 



AS vib ($) = S vib (<f>,T) - S mb ^o,T) =5S($- $ ), 



(11) 



and correspondingly, 



F vib (3>, T) = F OTb ($ , T) - T5S($ - $ c 



(12) 



where $o is a reference basin energy. The latter expression follows since the internal energy 
Uyi b = 3NksT for all basins. 

In addition to the basin free energy, the partition function in Eq.(|l|) requires knowledge 
of the configurational entropy density S c . Various recent studies have explored methods for 



estimating S c from computer simulations p7| , pO] , pT| , p8| , p!7| , p!8| . It has been observed that the 
distribution f2($) is well described by a Gaussian p?| , |29| , p!7| (equivalently, «S C ($) an inverted 
parabola). Although the arguments |27|,[25| m ay not apply to low energy minima, a Gaussian 
form for f2(<&) allows for a straightforward evaluation of the partition function Eq.(|l|), and 
whose validity has been tested in the range of temperatures where simulations are performed 
TT| , |2"2"|] . The configurational entropy density is written as 

5 C ($) ($ - $ ) 2 



Nk, 



a: 



(7" 



(13) 



where a is the height of the parabola and determines the total number of configurational 
states, i. e. energy minima (the total number is proportional to exp(aiV) ), <3> and o 2 
respectively define the mean and the variance of the distribution. The parameters a, $o 



and a have been estimated from simulation data [T7J. With the above form for «5> c ($) 
and Eq.(|T2|) for the vibrational free energy, the partition function can be evaluated, from 
which the following temperature dependence of the configurational entropy, the ideal glass 
transition temperature Tk (defined by S c (Tk) = and the inherent structure energies are 
obtained: 



< *>^ = *°"-mhf> < 14) 

where if/' = $„ + ^g, 

Ts, m = K&m err, - 1); m) = + ^) (i + T f) - gf B , 

(15) 

and 

T K = a(2Nk B Va + (r5S)- 1 . (16) 

These equations constitute relations that express quantities relevant to the thermodynamics 
of glass forming liquids, the configurational entropy and the ideal glass transition tempera- 
ture, in terms of paramters that describe the "energy landscape" of the liquid, namely the 
distribution of local energy minima, and the topography of individual minima in the form 
of vibrational frequencies. In particular, the expressions for TS C shows that the fragility of 
the liquid can be expressed in terms of parameters that quantify the "energy landscape" of 
the liquid. 



III. SIMULATION DETAILS AND RESULTS 



Many recent computer simulation studies of dynamics and thermodynamics aspects of 
glass forming liquids have employed a binary mixture of atomic particles as a model system 
|30| , |3llp3| , |20|JlO| , p!6|| , as this system has been parametrized to prevent crystallization. Results 
presented here are from molecular dynamics simulations of 204 type A and 52 type B 
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particles. The particles interact via the Lennard- Jones (LJ) potential, with parameters 
£ab/caa = 1.5, e BB /e AA = 0.5, o AB ja AA = 0.8, and o BB /a AA = 0.88, and m B /m A = 1. 
The LJ potential is modified with a quadratic cutoff and shifting at r^f = 2.ha a p P2"| . All 
quantities are reported in reduced units, length in units of a AA , volume V in units of a AA 
(density p = N/V, where N is the number of particles, in units of a AA = p ), temperature 
in units of e AA /k B , energy in units of e AA and time in units of r m = [a\ A mj 'e AA ) l l 2 , where 
m = m A = m B is the mass of the particles. Argon units are used for the A type particles 
when it is desirable to state values in SI units, i. e. e AA = 119. 8K x k B , a AA = 0.3405 nm, 
m A = 6.6337 x 10~ 26 kg. Molecular dynamics simulations are performed over a wide range 
of temperatures at each density, with run lengths ranging from 1.3ns to 0.4/zs. Local energy 
minimizations are performed for 1000 (k B T/e AA < 1.) or 100 (k B T/e AA > 1.) configurations 
to obtain typical local energy minima or 'inherent structures' sampled by the liquid. The 
Hessian (matrix of second derivatives of the potential energy) evaluated at the minima are 
diagonalized to obtain the vibrational frequencies. 

As the binary LJ liquid has mostly been studied at low temperatures and high densities, 
the location of the liquid-gas critical point has not previously been estimated. A rough 
estimate is obtained here by performing simulations in the appropriate range of p and T 
of length 2.912ns. Isotherms close to the critical point are shown in Fig. 1, which result 
in an estimate of the liquid-gas critical point of k B T c /e AA ~ 1.1, p c / Po ~ 0.416. The glass 
transition locus was estimated in |16| using diffusivity data fitted to the VFT form (Eq.(||)), 



and by thermodynamic means via Eq. (|4]) (by the condition that S c vanishes at the glass 
transition), where the total entropy of the liquid was obtained by thermodynamic integration 
from the ideal gas reference state (details in [0). The liquid-gas spinodal locus was also 
obtained to study the relation between these two limits to the liquid state. The two loci 
are seen to intersect at a finite T leading to the prediction of an ideal glass-gas mechanical 
instability locus at low temperatures. Fig. 2 shows this "phase diagram" of the binary LJ 
liquid, along with the locus of zero pressure (which shows that at low enough density, the 
liquid is under tension at the simulated temperatures). The figure also shows the T, p points 

8 



where simulations have been performed (see [0 for further details). 



The method of determining basin free energy and entropy (Eq.s from vibrational 
frequencies relies on the assumption that the basins are harmonic. Some indirect tests of 
this assumption have been presented recently elsewhere p^| . A direct test |15[ would be to 
thermally excite the inherent structures to excitation temperatures T e , and examine whether 
the energy of excitation above the inherent structure energy obeys the harmonic expectation 
of $(T e ) — $ m j„ = ^NksT. The result of this test is shown in Fig. 3 for p = 1.2p , where 
inherent structures are selected from simulations at three different temperatures, and are 
subjected to thermal excitations at T e that span a range around the simulation temperature. 
In all cases, the deviation from the harmonic expectation is less than 4%. Surprisingly, the 
deviations are bigger at lower temperatures. The reason is most likely that very closely 
related minima with negligible barriers exist at all energies studied |[23|| , and transitions 
between them occur even for low T e . Nevertheless, the results in Fig. 3 show that in 
practice, the harmonic approximation is quite valid for the binary LJ liquid. 

As described in the previous section, the classical calculation of the basin entropies in- 
dicates that differences in basin entropy between two basins is independent of temperature, 
and remain finite at zero temperature. Since this is not true in real materials (where all 
entropies vanish at zero temperature), one must examine to what extent conclusions from 
such calculations apply to real systems. Thre temperature range relevant to studying glass 
forming liquids is clearly the range above the glass transition temperature. Further, the 
relevant quantity for the analysis above was the difference of basin entropies. Therefore, if 
the quantum mechanical basin entropy differences correspond to a substantial fraction of the 
classical basin entropy differences above the glass transition, we may conclude that the clas- 
sical approximation is reasonable. Fig. 4 shows classical and quantum mechanical entropies 
obtained for basins sampled at two different temperatures (one 'high', one 'low') at two den- 
sities. It is clear (panels (a) and (d)) that the classical expression for the entropy deviates 
strongly from the quantum analog at low T (indeed, it diverges to — oo as T — > 0) while 
approaching it at high T. Similarly the quantum mechanical entropy differences approach 



the constant classical value at high T but decrease towards zero at low T. Nevertheless, at 
the glass transition temperature T Kl the quantum mechanical difference is still about 67% 
of the classical result. Since Argon units are used for the A atoms, and most glass formers 
are made of considerably bulkier molecules than atomic Argon, we may conclude that the 
classical result must apply well in most cases of interest. 



IV. SUMMARY AND CONCLUSIONS 



The discussion and results presented here indicate that the inherent structure approach 
is very useful in addressing interesting questions concerning the thermodynamics of glass 
forming liquids. It has been shown that the thermodynamic estimate of the glass transition 
locus agrees quite well with the one obtained from dynamic data [jlB] . Fruther, quantities 
such as the fragility of a liquid are reasonably well estimated by thermodynamic means, 
leading to an expression of fragility in terms of parameters that quantify the energy land- 
scape of the liquid. Additional data are presented here supporting the evaluation of basin 
entropies in a harmonic approximation and the classical calculations typically employed in 
such studies. 
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FIG. 1. Isotherms at temperatures close to the liquid-gas critical point, leading to an estiamate 
of the critical temperature and density, ksTc/tAA ~ 1-1, Pel 'po ~ 0.416. 



11 



1.05 



0.85 







(0 
a- 



Q. 

E 





0.65 



0.45 



0.25 



0.05 



\ N 



o o 



\ \ 
\ 



\ 



\ \ 
\ 



\ 



\ 



\ 



\ 



\ \ 
\ 



\ 



■ T from REMC 

s -i 

• T s from K T 

♦ T from fit to P 



\ ^ \ 

\ \ 

\ <m 

\ \ 

s ■- - ■ \ \ 

- T s Therm. Int. \ 
♦ T from Diffusivity ^ 

T K (Tl) \ \ 

*P* = 1.08p \ \ 

- T s (p-.08p )Therm. Int. \ 
Pressure = 0.0 V 

j x^L i \ lA 




0.4 



0.6 



0.8 1 

Density, p/p 



1.2 



o 



1.4 



FIG. 2. Liquid-gas spinodal obtained from (a) restricted ensemble Monte Carlo simulations 
(REMC), (b) extrapolations of isothermal compressibility kx, (c) polynomial fits to isotherms 
obtained in simulations, and (d) the empirical free energy obtained from thermodynamic integration 
('T s Therm. Int.'). The same curve is also shown shifted in p by 0.08/>o ( l T s (p— .08) Therm. Int.'). 
The glass transition locus obtained from (e) VFT fits to diffusivity data, and (d) extrapolation 
of configurational entropy to zero ('Tjg Therm. Int.'). (*) marks the density p* where inherent 
structure pressure is a minimum. The dot-dashed line represents the locus of zero pressure. Open 
circles represent T, p values where simulations have been performed at high densities. 
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FIG. 3. Test of the harmonic structure of minima. Energies shown are obtained by thermally 



exciting inherent structures from runs at different temperatures (shown in the legend), at 'excitation 
temperature's T e indicated on the x-axis and measuring the average potential energy. 
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FIG. 4. Basin entropies (panel (a) and (c)) and basin entropy differences (panels (b) and (d)) at 
two densities of the liquid (p = l.lpo for (a),(b), and l.2po for (c), (d)). At each density, two basins 
sampled at a low ("T") and high ("T n ") temperature are selected, for which the basin entropy 
vs. T is calculated both classically and quantum mechanically. For l.lpo, ksTi/eAA = 0.295 
ksTu/tAA = 0.6, and for 1.2po, ksTi/eAA = 0.47 ksTi/eAA = 0.8. Data shown indicates that in 
the quantum mechanical case, the entropy differences to indeed vanish as T — > 0, but at Tk, the 
quantum mechanical entropy difference is roughly 67% of the classical result. 
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